Non-equilibrium mechanics and dynamics of motor-activated gels 
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The mechanics of cells is strongly affected by molecular motors that generate forces in the cellular 
cytoskeleton. We develop a model for cytoskeletal networks driven out of equilibrium by molecular 
motors exerting transient contractile stresses. Using this model we show how motor activity can 
dramatically increase the network's bulk elastic moduli. We also show how motor binding kinetics 
naturally leads to enhanced low-frequency stress fluctuations that result in non-equilibrium diffusive 
motion within an elastic network, as seen in recent in vitro and in vivo experiments. 

PACS numbers: 87.16.Ka, 87.15.La, 62.20.Dc 



The mechanics of living cells are largely governed by 
the cytoskeleton, a complex network of filamentous pro- 
tein aggregates and various specialized proteins and en- 
zymes that couple the filaments together and gener- 
ate forces [![. As materials, in vitro networks of cy- 
toskeletal filaments have been shown to have unusual me- 
chanical properties, including a highly non-linear elastic 
response 0, y, |3, [H, Hi and negative normal stresses Q. 
Cytoskeletal networks in vivo, however, are far from 
equilibrium materials, due in large part to molecu- 
lar motors that exert internal forces within the net- 
works. This presents a challenge for quantitative sta- 
tistical/thermodynamic modeling. Recent studies of in 
vitro networks that include molecular motors have shown 
nearly a 100-fold stiffening of the networks due to mo- 
tor activity, as well as pronounced low-frequency, non- 
equilibrium fluctuations |8y. Here, we develop a model for 
such active gels that can explain both the strong stiff- 
ening of networks with motor activity, as well as the 
large non-equilibrium fluctuations at low frequencies. We 
also show how motor (un)binding kinetics naturally leads 
to a very simple and general form of stress fluctuations 
and diffusive-like motion, which are consistent with ob- 
served non-equilibrium dynamics in living cells 0, [Io| . 
This model can form the basis for quantitative design 
principles for creating synthetic polymeric materials with 
tunable elastic properties and muscle-like activation. 

Active solutions consisting of polymers and motors mo- 
tors constitute a strikingly new kind of material that can 
actively change/adapt its macroscopic mechanical prop- 
erties due to small-scale motor activit y th at drives rela- 
tive sliding of polymers past each otherjlll [l2l . [13 , lb! . [l5j . 
In permanently cross-linked networks, however, such 
motor activity can produce tensile stresses [8]. This 
muscle-like contraction is sketched in Fig. [TJ\. It is well 
known that single semi-flexible polymers stiffen under 
extension [III, and that this can result in macroscopic 
stiffening of networks under external strain 0, 0, [TtJ • This 



effect can also account for the observed dramatic stiff- 
ening of active networks d, [l8j]. Assuming an average 
state of tension in the network strands due to motor ac- 
tivity, we can calculate the expected degree of network 
stiffening as follows. The tension r in a single filament 
is calculated as a function of longitudinal extension I 
as in Ref. [171 ]. from which an effective spring constant 
K = dr/dt is calculated. In the nonlinear regime, this 
increases as K cx r 3 / 2 Q • The network modulus is given 
by G = j^p£ c K, where p is the density (length per vol- 
ume) of polymer, and £ c is the distance between cross- 
links [HI, [2(|. The predicted stiffening is shown in Fig. 
[1)3, where the filament tension has been normalized by 
the characteristic tension tq = kTir 2 i p / 'l\ required to 
pull out the fluctuations on a filament of length l c in 
the network. Here, l v is the persistence length. For a 
network of actin filaments, such as in Mizuno et al.Q, 
where i p = 17/xm and l c ~ 3/im, this characteristic av- 
erage tension is of order O.lpN, meaning that a tension 
of just a few pN, which is easily reached by myosin mo- 
tors, can lead to the observed 100-fold stiffening of active 
networks. 

The quasistatic picture sketched in Fig. [TK shows a 
motor (myosin minifilament) generating a pair of equal 
and opposite forces =pf applied at points f± = r ± a/2, 
separated by a. We expect a to be a few microns in 
an in vitro network. Since actin filaments are not able 
to support compressive loads over this distance, the re- 
sulting force dipole is contractile: the points are pulled 
together by a sort of muscle-like activity. While indi- 
vidual myosin motors are non-processive and are inca- 
pable of persistent, directed motion, they self-assemble 
into minifilamcnts, which arc processive. These minifil- 
aments still have a finite duty ratio. When they unbind 
the tension is instantaneously released, as sketched in the 
inset of Fig. [2|8| ]. Such a step-like force f(t) corresponds 
to a power spectrum of force fluctuations that varies as 
lu~ 2 , proportional to the square Fourier transform of /. 
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FIG. 1: (color online) (A) Schematic diagram of contractile 
motor activity in a network. A myosin minifilament (blue) 
slides two network filaments (red) past each other, generating 
an equal and opposite pair of forces (green arrows). (B) Plot 
of the predicted relative stiffening of a semi-flexible network 
as a function of (normalized) motor-induced tension. The 
inset shows the nonlinear force-extension relation of a single 
semi-flexible filament 0, 0, Q3 • 



As we show, this physical picture of step-like contrac- 
tile forces naturally leads to non-equilibrium fluctuations 
that dominate only at low frequencies, as sketched in Fig. 
[U Surprisingly, this generates motion that appears to be 
diffusive: (\x(t)—x(0)\ 2 ) ~ Dt, but occurring in an elastic 
material. The effective diffusion constant D is controlled 
by motor activity and not temperature. Using well- 
established viscoelastic properties of cross-linked F-actin 
networks [TH |2(j, we find distinct regimes of both ther- 
mal and athermal (motor-induced) fluctuations sketched 
in Fig. [21 which are consistent with the observations both 
in vivo ^[3] and in vitro^. 

To model the active gel we use a continuum description 
for a viscoelastic homogeneous and isotropic medium, but 
in which the motor activity couples to this medium as il- 
lustrated in Fig. [T]A_. For in vitro networks such as in 
Ref. [1], the distance between cross-links, and thus a is 
expected to be of order 3-fO^m. On this scale, we can 
model the action of a motor as the introduction of a pair 
for equal and opposite applied forces in the (visco-) elastic 
continuum. The resulting displacement field u, at posi- 
tion fo of the network we describe by a linear response 
function depending on position and frequency as 

Ui (rb, lu) = [otij (fo -f-,u)- a i:j (f Q - r + , oj)] f 3 (w), 
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FIG. 2: The displacement power spectral density (PSD) in an 
active gel. Here, frequency is measured in terms of £1 = uT/B. 
The thermal PSD (dashed line) shows a plateau at low fre- 
quencies. Thus, the active component of the PSD dominates 
at low frequencies, while the thermal PSD is expected to 
dominate at high frequencies. (Inset) Schematic of the time- 
dependent force due to molecular motor activity. 



using the fact that the motor-generated forces =f/ are 
equal and opposite. Stability also requires that / and a 
be parallel. The response function to a point force ay- 
can be written in terms of a\\ and aj_, where (f) — 
hfja\\ (r) + (Sij - fifj) a± (r). 

We calculate these two response components within a 
two-fluid approximation, in which the cytoskclctal fila- 
ments are treated as a porous elastic network immersed 
in a viscous solvent [U [H, H3, 0|- Here, the network 
displacement u and solvent velocity v satisfy the coupled 
equations 



= /iV 2 u +(p + A)V(V -u)+T 



= rjW v — VP - r 



du 



+ fs 



(2) 



(3) 



where \i and A are Lame coefficients, r\ is the solvent vis- 
cosity, and the forces f n}S represent the forces on the net- 
work and solvent, respectively. Given a mcshwork with 
a pore size £, the coupling T is expected to be of order 
7//£ 2 . These are solved for the response of the combined 
system to an applied point force. The resulting response 
functions are given by 
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e -2 ^-' _ x±(x). Here, G is the shear modulus and 
B = 2{ \-2a G ^ S ^ ne longitudinal modulus, where a is 
the Poisson ratio, and = ujY/B. This coupling can 
be understood in terms of the solvent flow through the 



3 



highly porous gel: rapid solvent flow through the fila- 
ment mesh gives rise to large shear stresses, effectively 
dragging the network with the solvent. This drag pre- 
vents the large-scale relative motion of the network and 
solvent beyond a range of order Or 1 / 2 . On larger length 
scales r or at higher frequencies u>, the drag effectively in- 
hibits the relative motion of solvent and network so that 
for r^/fi 3> 1, the combined network and solvent act as 
a single incompressible material [23l |24| . and X\\,± both 
vanish (Fig. EJ^.). Here, the response of the medium is 
purely transverse (the displacement vector field is diver- 
genceless) and is given by the generalized Oseen tensor, 
given by leading terms in square brackets above [24]. The 
corresponding volume-preserving flow response of an in- 
compressible gel when subject to a symmetric pair of 
point forces is shown in Fig. HJA. 

In this incompressible case, the displacement field u(uj) 
of the network resulting from motor activity varies with 
an overall frequency dependence proportional to the ratio 
of the force f(oj) to the shear modulus G(w), according to 
Eqs. (1-3). Thus, we find for the model illustrated in Fig. 
m that (\u(cj)\ 2 ) oc {\f{uj)\ 2 )/\G(uj)\ 2 oc \luG\~ 2 . Cross- 




FIG. 3: (A) Graphs of the spatial dependence of the longitu- 
dinal parts of the parallel (||) and perpendicular (_L) response 
functions (Eqs. (2) and (3)). The effect of compression of 
the network on the response functions can be reduced to a 
universal form when p lotted against the dimensionless quan- 
tity r\/Vi — r^JuoT / B, demonstrating the diffusive nature of 
the propagation of the network density mode. (B) The effect 
of network compression can be isolated in experimental data 
by examining the difference in the parallel and perpendicular 
response functions given in Eq. (4). Here we plot the pre- 
dicted form of the real (Re) and imaginary (Im) parts of that 
difference vs. the dimensionless variable ry/Ti. 



linked biopolymer networks typically exhibit a constant 
or weakly frequency-dependent elastic regime as a func- 
tion of frequency. Here, we expect to see (\u(co)\ 2 ) oc w~ 2 , 
which is consistent with recent displacement fluctuations 
observed in cells|9|], and which corresponds to diffusive 
motion. At higher frequencies, such networks typically 
exhibit a power-law increase in the shear modulus with 
frequency |19l. [20L l23j ], in which G oc w 3 / 4 . In this fre- 
quency regime stain fluctuations in the active gel take 



the form (\u(lu) 



OC OJ 



-7/2 



as shown in Fig. [2] For com- 



parison, the equilibrium thermal fluctuations for such a 
network are shown as the dashed line. At low frequen- 
cies the motor-driven fluctuations will dominate over the 
ever-present thermal fluctuations, consistent with the re- 
sults of both Lau et al. § and Mizuno et al. [|[ . 

Since biopolymer and cytoskeletal networks are gener- 
ically porous with pore sizes of order 1/zm, they can 
deform compressibly. This density mode, however, is 
strongly suppressed by drag at high enough frequencies. 
The loss of the density mode at high frequencies is illus- 
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FIG. 4: (A) The displacement vector field of an incompress- 
ible network shown in a plane passing through the two force 
centers for a contractile motor acting at the origin. The forces 
are applied symmetrically at points (±3/4, 0) and are each di- 
rected towards the origin. (B) The network displacement field 
for the compression mode shown in the limit of low frequency 
or weak hydrodynamic coupling (r — > 0). Again the forces 
are applied symmetrically at points (±3/4,0) and are each 
directed towards the origin. The resulting displacement field 
induces network density variations in the material. 
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trated in Fig. [3j\, where the effects of finite compressibil- 
ity, represented by X\\,X> vanish at high frequency. Al- 
though the basic physics of these effects have been dis- 
cussed before for both flexible polymer systems [2ll [22| 
and semi-flexible biopolymer systems [23j, |24|, there has 
been no direct experimental observation of these com- 
pressibility effects in porous biopolymer systems. 

We can isolate the effects of the network compressibil- 
ity by the examining the combination 

! [xx (rVn) - |x|| (»Vn) 
a±(r, U ) - -on(r, W ) = , 

(6) 

which is plotted in Fig.[3p. This measurable combination 
of response functions strictly vanishes in the incompress- 
ible limit. This, along with the specific combined r and lj 
dependence, may permit the first direct measurement of 
compressibility effects that are expected to be character- 
istic of biopolymer/cytoskeletal networks. Furthermore, 
the flow/displacement field corresponding to this com- 
pressible mode (shown in Fig. 2J3 in the limit V — > 0) 
strongly differ from the case of an incompressible sys- 
tem (Fig. SK). Here, the longitudinal (irrotational dis- 
placement field) contributions to the response function 
are an = and — l/(8ivrB). The difference in 
spatial structure of these strain fields may also be used 
to experimentally identify the effects of compression. 

To consider the effect of multiple contractile events 
within the medium, we can represent the resulting dis- 
placement field at the origin itj at by a sum 

Ui = ^2 Aa u o) %7> ( 7 ) 

where Acty (r, a) = onj (f— a/2) — aij (f+ a/2) is the 
response to a contractile force pair. We suppress the fre- 
quency dependence. This sum represents the combined 
effect of temporally uncorrelated contractile events oc- 
curring homogeneously throughout the medium. This 



assumption remains valid provided that the events rarely 
occur with a separation of order a ~ t c during the typi- 
cal processivity time to- Such a sum or average has been 
performed in calculating the PSD in Fig. [2] for the case 
of an incompressible network. In this case the scaling 
described above is a good approximation. 

This model shows how motor activity within a semi- 
flexible gel, together with the well-established non-linear 
response of such networks leads to a strong stiffening 
of the network, and that this stiffening increases more 
than linearly with the motor force. This can account for 
the recently observed nearly 100-fold network stiffening 
with motor forces of order 1-10 pN[8(. Furthermore, the 
(un)binding kinetics of the motors naturally leads to a 
specific characteristic time dependence of the force fluc- 
tuations in active gels. Given a finite processivity time 
to over which minifilaments remain bound and generate 
force, the unbinding results in 1 /uj 2 force fluctuations for 
frequencies uj > 1/to- This spectrum is a direct result of 
the expected sharp time dependence of motor unbinding, 
and is insensitive to slow variations of force during motor 
motion. For frequencies u) > 1/io, the divergence of the 
force spectrum will be suppressed. Our model, is for un- 
correlated motor activity, in that the total fluctuations 
can be represented as a sum of independent fluctuations 
due to individual motor force generation and unbinding. 
At sufficiently high motor densities, one might expect co- 
operativity of motor activity, whose consequences can be 
studied in extension of the present model. 
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